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Mitochondrial  DNA  analysis  is  very  useful  for  the  interpretation  of  the  history  of  human  migration  and  to 
estimate  the  frequency  of  a  haplotype  in  the  forensic  context.  From  a  human  settlement  perspective,  La 
Paz  area  is  greatly  interesting  since  the  first  planned  city  of  the  region  is  located  there.  Samples  from  110 
individuals  from  La  Paz  were  studied  analysing  the  polymorphisms  in  the  D-loop,  hypervariable  region  I 
(HVI)  and  hypervariable  region  II  (HVII)  in  order  to  verify  the  genetic  diversity.  The  aim  of  this  study  was 
to  start  the  creation  of  a  population  database  in  order  to  obtain  the  genetic  interpopulation  variability 
and  classify  haplotypes  into  characteristic  haplogroups  of  South  America.  A  total  of  97  different  haplo- 
types  were  identified,  90  being  unique,  expressed  by  122  polymorphic  nucleotide  positions.  Nucleotide 
and  sequence  diversity  were  estimated  to  be  0.015  ±  0.0075  and  0.996,  respectively.  Haplogroup 
distribution  in  the  samples  was  57.27 %  B4, 19.09%  Cl,  10.00%  A2,  3.64%  Dl,  2.73%  D4h3, 1.82%  H,  and  0.91% 
for  each  of  the  haplogroups  A4,  B4cla,  CZ,  D4J,  M7a  and  M8/N9b.  The  rate  of  length  heteroplasmy  was 
36.36%  in  HVI  and  52.73%  in  HVII.  Phylogenetic  analysis  reveals  proximity  to  the  Korean,  Chilean 
aboriginal,  Japanese  and  Australian  populations.  The  estimated  genetic  variability  of  the  studied  pop¬ 
ulation  was  high,  suggesting  an  early  settlement. 

©  2010  Elsevier  Ltd  and  Faculty  of  Forensic  and  Legal  Medicine.  All  rights  reserved. 


1.  Introduction 

Bolivia  is  a  South  American  country  in  the  heart  of  the  Andes. 
The  history  of  the  territory  follows  the  cultures  of  the  first  empire 
within  the  domains  of  the  first  planned  city  of  the  region,  Tiwa- 
naku,  sited  in  the  South  coast  of  the  Lake  Titicaca.  The  last  continent 
to  be  populated  by  humans  was  the  Americas  and  several  questions 
remain  open  regarding  this  event.  Due  to  its  antiquity,  the  Andean 
region  can  reveal  some  of  the  human  migration  and  colonization 
history  of  the  New  World  since  most  of  the  Bolivian  population 
descends  from  Aymaras  and  Quechuas,  two  cultural  groups  from 
the  Tiwanaku  and  Inka  Empire,  which  persisted  labouring  the  fields 
of  the  Andean  highlands  around  Lake  Titicaca  until  the  present 
days.1 

Mitochondrial  DNA  (mtDNA)  is  an  excellent  tool  for  population 
studies  and  forensic  genetics  due  to  its  high  copy  number  per  cell, 
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maternal  inheritance,2  high  mutation  rate  of  evolution  that  provide 
highly  polymorphism  to  the  sequences  and  absence  of  recombi¬ 
nation.3  The  discrimination  power  of  mtDNA  is  due  to  the  poly¬ 
morphic  nature  of  the  hypervariable  regions  (HVI,  HVII,  and  HVIII) 
which  are  located  in  the  control  region.  The  mitochondrial  genome 
has  become  the  most  widely  used  genetic  marker  in  human 
evolution  studies.  Sequences  with  shared  mutations  that  tend  to 
reveal  a  regional  specificity  define  haplogroups.  These  haplogroups 
are  important  to  clarify  the  history  and  demographic  past  of 
a  population  since  they  can  reflect  a  phylogenetic  relationship 
between  populations.4  The  first  mtDNA  haplogroups  discovered  in 
Native  Americans  were  determined  by  Torroni  et  al.,  1993.5  Poste¬ 
rior  studies  of  mtDNA  variation  concluded  that  Native  American 
populations  have  five  distinct  major  haplogroups,  all  of  Asian  origin 
(A,  B,  C,  D  and  X).5,6  For  forensic  purposes,  the  discrimination  power 
between  individuals  based  on  mtDNA  is  lower  than  traditional 
nuclear  DNA  markers  because  of  its  maternal  inheritance;  however, 
mtDNA  analysis  is  particularly  useful  when  there  are  degraded  or 
low  level  DNA  samples.  To  solve  cases  of  missing  persons  and  to 
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Table  1 

HVI  and  HVII  mtDNA  sequences  of  110  individuals  from  La  Paz.  N  —  number  of  occurrence  of  the  haplotype;  Est.  Hg  —  “estimated  haplogroup”. 


HVI 


HVII 


LPAZ001 

16183C; 

LPAZ002 

16223T; 

LPAZ003 

16183C; 

LPAZ004 

16093C; 

LPAZ005 

16183C; 

LPAZ006 

16183C; 

LPAZ007 

16183C; 

LPAZ008 

16183C; 

LPAZ009 

16223T; 

LPAZ010 

16183C; 

LPAZ01 1 

16183C; 

LPAZ012 

16183C; 

LPAZ013 

16182C; 

LPAZ014 

16183C; 

LPAZ015 

16168T; 

LPAZ016 

16223T; 

LPAZ017 

16168T; 

LPAZ018 

16183C; 

LPAZ019 

16182C; 

LPAZ020 

16183C; 

LPAZ021 

16183C; 

LPAZ022 

16183C; 

LPAZ023 

16183C; 

LPAZ024 

16183C; 

LPAZ025 

16183C; 

LPAZ026 

16183C; 

LPAZ027 

16183C; 

LPAZ028 

16129A; 

LPAZ029 

16183C; 

LPAZ030 

16183C; 

LPAZ031 

16183C; 

LPAZ032 

16183C; 

LPAZ033 

16183C; 

LPAZ034 

16223T; 

LPAZ035 

16183C; 

LPAZ036 

16183C; 

LPAZ037 

16223T; 

LPAZ038 

16105K; 

LPAZ039 

16183C; 

LPAZ040 

16183C; 

LPAZ041 

16066G; 

LPAZ042 

16114G; 

LPAZ043 

16183C; 

LPAZ044 

16051G; 

LPAZ045 

16183C; 

LPAZ046 

16192T; 

LPAZ047 

16129A; 

LPAZ048 

16111T; 

LPAZ049 

16183C; 

LPAZ050 

16182C; 

LPAZ051 

16183C; 

LPAZ052 

16223T; 

LPAZ053 

16176T; 

LPAZ054 

16092C; 

LPAZ055 

16051G; 

LPAZ056 

16092C; 

LPAZ057 

16182C; 

LPAZ058 

16223T; 

LPAZ059 

16223T; 

LPAZ060 

16223T; 

LPAZ061 

16223T; 

LPAZ062 

16172C; 

LPAZ063 

16111T; 

LPAZ064 

16223T; 

LPAZ065 

16223T; 

LPAZ066 

16037G; 

LPAZ067 

16183C; 

LPAZ068 

16093C; 

16298C; 

LPAZ069 

16111T; 

LPAZ070 

16183C; 

LPAZ071 

16092C; 

LPAZ072 

16140C; 

16189C;  16217C3 
16292T;  16325C;  16362C 
16189C;  16217C3 
16183C;  16189C;  16217C3 
16188T;  16189C;  16217C 
16188T;  16189C;  16217C 
16188T;  16189C;  16217C 
16189C;  16217C3 
16325C;  16362C 
16189C;  16217C3 

16189C; 16217C;  16270T;  16278T3 
16188T;  16189C;  16217C; 16362C 
16183C;  16189C;  16217C3 
16189C;  16217C3 

16183C;  16189C;  16217C; 16270T;  16278T3 
16241G;  16304C;  16325C;  16362C 
16182C;  16183C;  16189C;  16217C;  16218T3 
16188T;  16189C;  16217C 
16183C;  16189C;  16217C;  16362C3 
16189C;  16217C;  16293C3 
16189C;  16217C3 
16188T;  16189C;  16217C;  16354T 
16189C;  16217C3 
16188T;  16189C;  16217C; 16266T 
16188T;  16189C;  16217C;  16354T;  16362C 
16189C;  16217C3 

16189C;  16217C;  16261T;  16319A3 
16183C;  16189C;  16217C3 
16189C;  16217C;  16289G3 
16189C;  16217C;  16356C;  16362C3 
16189C;  16217C3 
16189C;  16217C;  16359C3 
16188T;  16189C;  16217C; 16354T 
16255A;  16301T;  16342C;  16362C 
16189C;  16217C;  16360T3 
16189C;  16217C;  16274A;  16319A; 16362C3 
16266T;  16298C;  16325C;  16327T 
16168T;  16183C;  16189C;  16217C;  16301A3 
16189C;  16217C3 

16188T;  16189C;  16217C;  16304C;  16362C 
16183d; 16186T;  16189C;  16217C 
16129A; 16222T;  16223T;  16325C;  16362C 
16189C;  16217C;  16289G3 
16182C;  16183C;  16189C;  16217C3 
16189C;  16217C;  16301A3 
16223T;  16298C;  16325C;  16327T 
16223T;  16298C;  16325C;  16327T 
16223T;  16290T;  16319A;  16362C 
16189C;  16217C;  16319A; 16360T3 
16183C;  16189C;  16217C;  16293G; 1631 1C3 
16189C;  16217C;  16362C3 
16298C;  16325C;  16327T 

16183C;  16189C;  16217C;  16240G;  16289G; 16362C3 

16145A;  16223T;  16298C;  16325C 

16129A; 16182C; 16183C; 16189C; 16217C; 16234A3 

16183C;  16189C;  16217C;  16362C3 

16183C;  16189C;  16217C;  16362C3 

16293G;  16325C;  16362C 

16290T;  16319A;  16362C 

16298C;  1631 1C;  16325C;  16327T 

16298C;  16325C;  16327T 

16223T;  16256T;  16298C;  16325C; 16327T 

16217C;  16223T;  16290T;  16319A;  16343T;  16362C 

16298C;  16325C;  16327T;  16345T 

16298C;  16325C;  16327T;  16381C 

16223T;  16298C;  16325C;  16327T 

16189C;  16217C;  16270T;  16278T3 

16183C;  16189C; 16223T; 

16325C;  16327T;  16381C3 

16217C; 16223T;  16290T;  16319A;  16343T;  16362C 
16188T;  16189C;  16217C 

1611 IT;  16223T;  16290T;  16319A;  16362C 
16223T;  16298C;  16311C;  16325C;  16327T 


73G;  207A;  263G;  315.1C 

73G;  263G;  315.1C 

73G;  103 A;  263G;  315.1C 

73G;  263G;  309.1C;  315.1C 

73G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  185A;  263G;  309.1C;  315.1C 

73G;  241G;  263G;  309.1C;  315.1C3 

73G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  143 A;  263G;  309.1C;  315.1C 

73G;  146C;  263G;  309.1C;  309.2C;  315.1C3 

73G;  263G;  309.1C;  315.1C3 

73G;  152C;  186T;  263G;  315.1C 

73G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  143 A;  146C;  215G;  263G;  315.1C 

73G;  263G;  315.1C 

73G;  263G;  309.1C;  315.1C 

73G;  263G;  315.1C 

73G;  151T;  186T;  204C;  263G;  315.1C 
73G;  263G;  309d;  315.1C; 

73G;  204C;  207A;  263G;  31 5.1  C3 

73G;  151T;  263G;  309.1C;  309.2C;  315.1C3 

73G;  186T;  263G;  309.1C;  315.1C 

73G;  143 A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  186T;  263G;  309.1C;  309.2C;  315.1C3 

73G;  186T;  263G;  309.1C;  315.1C3 

73G;  103 A;  143A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  75A;  263G;  315.1C;  340T3 

73G;  146C;  263G;  309.1C;  309.2C;  315.1C3 

73G;  143 A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  204C;  263G;  309.1C;  315.1C3 

73G;  207A;  257G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  143 A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  146C;  186T;  263G;  309.1C;  315.1C 

73G;  152C;  263G;  309.1C;  315.1C3 

73G;  143 A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  263G;  309.1C;  315.1C3 

73G;  249d;  263G;  290d;  291d;  315.1C 

73G;  204C;  263G;  309.1C;  315.1C3 

73G;  143 A;  146C;  215G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  152C;  186T;  263G;  309.1C;  315.1C 

73G;  143 A;  210G;  263G;  309.1C;  315.1C 

73G;  143 A;  263G;  309.1C;  315.1C 

73G;  143 A;  146C;  215G;  263G;  309d;  315.1C 

73G;  152C;  263G;  309.1C;  309.2C;  315.1C3 

73G;  146C;  197G;  215G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  146C;  153G;  263G;  309.1C;  309.2C;  315.1C3 

73G;  146C;  150T;  263G;  309.1  C;  309.2C;  315.1C3 

73G;  143 A;  263G;  309.1C;  309.2C;  315.1C3 

73G;  103 A;  146C;  151T;  263G;  309.1C;  309.2C;  315.1C3 

73G;  79T;  249d;  263G;  290d;  291d;  309.1C;  315.1C3 

73G;  204C;  263G;  309.1C;  315.1C3 

73G;  185A;  249d;  263G;  290d;  291d;  315.1C 

73G;  152C;  263G;  309.1C;  315.1C3 

73G;  143 A;  146C;  215G;  263G;  309.1C;  315.1C3 

73G;  146C;  215G;  263G;  309.1C;  309.2C;  315.1C3 

97A;  106d;  107d;  108d;  109d;  llOd;  llld;  263G; 315.1C 

73G;  146C;  153G;  182T;  194T;  195C;  235G;  263G;  315.1C 

73G;  150T;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  162Y;  249d;  263G;  290d;  291d;  309.1C;  309.2C;  315.1C3 

73G;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  146C;  235G;  260A;  263G;  315.1C 

73G;  146C;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  152C;  249d;  263G;  290d;  291d;  309.1C;  315.1C3 

73G;  222T;  249d;  263G;  290d;  291d;  309.1C;  315.1C 

73G;  94 A;  152C;  183G;  204C;  263G;  309.1C;  315.1C3 

73G;  249d;  263G;  290d;  291d;  315.1C 

73G;  146C;  153G;  235G;  260A;  263G;  315.1C 

73G; 106d; 107d; 108d; 109d;  llOd; llld;  186T;  263G; 

309.1C;  315.1C 

73G;  146C;  153G;  222T;  235G;  263G;  309.1C;  315.1C 
73G;  249d;  263G;  290d;  291d;  309.1C;  309.2C;  315.1C3 


N  Est.  Hg 
1  B4 

1  D1 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  D1 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  D1 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  B4 

1  D4h3 

1  B4 

1  B4 

1  Cl 

1  B4 

1  B4 

1  B4 

1  H 

1  D1 

1  B4 

1  B4 

1  B4 

1  Cl 

1  Cl 

1  A2 

1  B4 

1  B4cla 

1  B4 

1  Cl 

1  B4 

1  CZ 

1  B4 

1  B4 

1  B4 

1  D4j 

1  A4 

1  Cl 

1  Cl 

1  Cl 

1  A2 

1  Cl 

1  Cl 

1  Cl 

1  B4 

1  Cl 

1  A2 

1  B4 

1  A2 

1  Cl 
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Table  1  ( continued ) 


HV1 


HV1I 


N  Est.  Hg 


LPAZ073  16223T;  16298C;  16325C;  16327T 

LPAZ074  16223T;  16298C;  16325C;  16327T 

LPAZ075  16145A; 16223T;  16270T;  16298C;  16325C;  16327T 

LPAZ076  1611  IT;  16136C;  16223T;  16268A;  16290T;  16319A;  16362C 

LPAZ077  1611 IT;  16217C;  16223T;  16290T;  16319A;  16343T;  16362C 

LPAZ078  1611 IT;  16183G; 16223T;  16298C;  16325C;  16327T;  16368C 
LPAZ079  16075C;  16223T;  16241G; 16260T;  16294T;  16301T; 

16342C;  16362C 

LPAZ080  16181G;  16189C;  16223T;  16298C;  16325C;  16327T;  16344Ta 
LPAZ081  161 1 IT;  16217C;  16223T;  16290T;  16319A; 16362C 
LPAZ082  16051G;  16213A;  16223T;  16247G; 16298C; 16325C;  16327T 
LPAZ083  16075C;  16223T;  16241G; 16260T;  16294T;  16301T; 

16342C;  16362C 

LPAZ084  16041G;  16172C;  16173T;  16192T;  16223T; 

16266T;  16298C;  16325C;  16327T;  16346A 
LPAZ085  161 1 IT;  16183C;  16189C;  16217C;  16223T; 

16290T;  16319A;  16343T;  16362Ca 
LPAZ086  16189C;  16192T;  16223T;  16266T;  16298C; 

16316G;  16325C;  16327T 

LPAZ087  16189C;  16192T;  16223T;  16266T;  16298C;  16316G; 

16325C;  16327T 

LPAZ088  16183C; 16188T;  16189C;  16217C 

LPAZ089  16183C; 16188T;  16189C;  16217C 

LPAZ090  161 1 IT;  16223T;  16290T;  16319A; 16362C 

LPAZ091  16183C; 16188T;  16189C;  16217C 

LPAZ092  16181G;  16189C;  16217C 

LPAZ093  16209C;  16223T;  16311C 

LPAZ094  16182C;  16183C;  16189C;  16223T;  16298C;  16325C; 

16327T;  16344Ta 

LPAZ095  16183C; 16188T;  16189C;  16217C 

LPAZ096  16140C;  16183C;  16188T;  16189C;  16217C 

LPAZ097  16183C; 16188T;  16189C;  16217C 


73G;  146C;  195C;  249d;  263G;  290d;  291d;  309.1C;  315.1C;  340T  1  Cl 

73G;  125C;  127C;  249d;  263G;  290d;  291d;  309.1C;  1  Cl 

309.2C;  315. lCa 

73G;  214G;  249d;  263G;  290d;  291d;  309.1C;  315.1C  1  Cl 

73G;  146C;  153G;  185A;  235G;  263G;  309.1C;  315.1C3  1  A2 

73G;  146C;  152C;  153G;  235G;  260A;  263G;  315. 1C  1  A2 

73G;  152C;  249d;  263G;  290d;  291d;  309.1C;  315.1Ca  1  Cl 

73G;  94A;  152C;  195C;  263G;  309.1C;  315.1C  1  D4h3 

73G;  185A;  249d;  263G;  290d;  291d;  309.1C;  315.1Ca  1  Cl 

73G;  146C;  153G;  212C;  215G;  235G;  263G;  309.1C;  309.2C;  315.1Ca  1  A2 

73G;  146C;  195C;  249d;  263G;  290d;  291d;  309.1C;  315.1C  1  Cl 

73G;  94A;  125C;  127C;  152C;  195C;  263G;  309.1C;  315.1C  1  D4h3 


73G;  204C;  249d;  263G;  290d;  291d;  315.1Ca 

73G;  146C;  153G;  235G;  260A;  263G;  309.1C;  309.2C;  315.1Ca 


1  Cl 
1  A2 


73G;  143 A;  150T;  249d;  263G;  290d;  291d;  309.1C;  315. 1C; 
325T;  338T 

73G;  143 A;  150T;  201G;  249d;  263G;  290d;  291d;  309.1C; 

315. 1C;  325T 

73G;  186T;  263G;  315.1C 

73G;  186T;  263G;  309.1C;  315.1C 

73G;  146C;  153G;  235G;  263G;  309.1C;  315.1C3 

73G;  263G;  309.1C;  315.1C 

73G;  185A;  249d;  263G;  290d;  291d;  309.1C;  315.1Ca 
73G;  150T;  152C;  189G;  195C;  200G;  263G;  309.1C;  315.1Ca 
73G;  263G;  309.1C;  309.2C;  315.1Ca 


1  Cl 

1  Cl 

2  B4 

2  B4 

3  A2 
5  B4 
1  H 

1  M7a 
1  M8/N9b 


263G;  309.1C;  315.1Ca  2  B4 

263G;  315. 1C  2  B4 

263G;  309.1C;  309.2C;  315.1Ca  4  B4 


a  Means  length  heteroplasmy. 


ascertain  the  identity  of  victims  of  mass  disasters  or  terrorism,  it  is 
important  to  have  a  high  quality  population  database  at  disposal 
that  allows  to  estimate  the  frequency  of  a  questioned  haplotype.7  At 
the  moment,  however,  mitochondrial  DNA  databases  are  limited  in 
representative  population  data  of  South  American  origin.  This 
paper  presents  a  contribution  to  the  worldwide  collection  of 
mtDNA  sequences  set  on  HVI  and  HVU  haplotypes  from  La  Paz  for 
forensic  and  phylogenetic  purposes. 

The  multiple  mtDNA  copies  in  one  individual  are  in  general 
identical,  a  condition  known  as  homoplasmy.  But  homoplasmy  is 
not  universal  and  the  co-existence  of  different  mtDNA  sequences  in 
one  individual  is  possible,  a  condition  identified  as  heteroplasmy. 
Heteroplasmy  is  more  frequently  detected  in  the  control  region  as 
length  variants,  named  length  heteroplasmy.  Point  heteroplasmy  is 
the  product  of  a  nucleotide  change  at  a  certain  position  that  some 
sequences  have  and  others  do  not.  It  has  been  determined  that,  at 
present,  the  rate  of  heteroplasmy  in  one  population  is  14%.8 

This  paper  reports  the  mtDNA  of  one  population  from  La  Paz 
(Bolivia)  composed  of  haplotypes  from  110  native  individuals,  the 
distribution  of  the  haplogroups  in  the  population  and  a  phyloge¬ 
netic  study  in  relation  to  other  populations. 

The  aims  of  this  work  were  to  (i)  start  the  creation  of  a  mtDNA 
database  of  the  La  Paz  population,  (ii)  obtain  the  mtDNA  variability, 
(iii)  classify  the  haplotypes  into  haplogroups  in  order  to  discover 
the  genetic  structure  of  La  Paz  population  through  mtDNA  and  (iv) 
infer  the  evolutionary  processes  that  have  formed  the  history  of  the 
La  Paz  population. 

2.  Material  and  methods 

Blood  samples  were  obtained  respecting  the  informed  consent 
established  in  Bolivia  from  110  healthy  and  unrelated  individuals 


with  ages  between  sixteen  and  seventy  two  years  old.  All  of  the 
individuals  are  from  the  La  Paz  department  of  Bolivia.  DNA  was 
extracted  using  the  Chelex®100  method.9  The  hypervariable  regions 
of  the  mtDNA  non-coding  control  region  were  amplified  using  the 
primers  L15997  and  H16401  for  HVI  and  the  primers  L408  and  H048 
for  HVI110  and  thermocycling  conditions  were  performed  in 
a  GeneAmp®  PCR  system  2700  (Applied  Biosystems,  Foster  City,  CA). 
PCR  amplification  was  carried  out  using  IX  Buffer  PE  11,  MgCl2  at 
2  mM;  dNTPs  at  0.2  mM  each,  primers  at  0.2  pM,  5  U/pl  of  AmpliTaq 
Gold®  DNA  Polymerase  (Applied  Biosystems,  Foster  City,  CA)  and 
2  ng  of  genomic  DNA  in  a  25  pL  final  reaction  volume.  Thermocycling 
conditions  was  as  follows:  at  95  °C  for  6  min,  followed  by  32  cycles  of 
10  s  at  95  °C,  30  s  at  60  °C,  30  s  at  72  °C,  and  finally  at  15  °C  for  1  min. 
Both  light  and  heavy  chains  of  each  hypervariable  region  were 
sequenced  using  the  sequencing  chemistry  ABI  Prism®  dRhodamine 
Terminator  Cycle  Sequencing  Ready  Reaction  Kit  (Applied  Bio¬ 
systems,  Foster  City,  CA).  The  sequencing  cycle  was  performed  in 
a  GeneAmp®  PCR  system  2700  (Applied  Biosystems,  Foster  City,  CA) 
in  a  25  cycles:  at  96  °C  for  10  s,  56  °C  for  5  s,  and  72  °C  for  4  min.  The 
nucleotide  sequence  of  HVI  from  position  16024  until  position 
16390,  and  nucleotide  sequence  of  HV1I  from  position  73  until 
position  340  were  analysed  in  the  ABI  PRISM®  sequencers  3100  - 
Avant  Genetic  Analyser  (with  the  DNA  Sequencing  Analysis 


Table  2 

Diversity  parameters  for  mtDNA  sequences  of  110  individuals  from  La  Paz. 


MtDNA  control  region 

Nucleotide  diversity  (p) 

Sequence  diversity  (h) 

HVI 

0.011 

0.948 

HV2 

0.015 

0.986 

HVI +  HV2 

0.015  ±  0.0075 

0.996 
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■  Haplogroup  B4  -  57.27% 


■  Haplogroup  Cl  - 19.09% 


■  Haplogroup  A2  - 10.00% 


■  Haplogroup  D1  -  3.64% 


■  Haplogroup  D4h3  -  2.73% 


■  Haplogroup  H  - 1.82% 


■  Other  Haplogroups  -  5.46% 


Fig.  1.  Haplogroup  distribution  of  mtDNA  in  La  Paz,  Bolivia. 


Software™  v.3.7  and  the  SeqScape®  Software  v.2.0)  and  in  the  3130  — 
Genetic  Analyser  (with  the  AB1  DNA  Sequencing  Analysis  Software 
v.5.2  and  the  SeqScape®  Software  v.2.5).  The  obtained  haplotypes 
were  compared  with  the  Cambridge  Reference  Sequence  (CRS)11,12 
and  the  differences  found  in  each  sequence  regarding  the  CRS 
were  typed  following  the  nomenclature  recommendations  of  the 
1UPAC.13  The  establishment  of  the  haplogroups  was  made  consid¬ 
ering  the  nucleotide  substitutions  on  non-coding  hypervariable 
regions  1  and  11  (HV1  +  HV1I)  on  a  Web-based  tool  for  management 
and  quality  analysis  of  mitochondrial  DNA14  based  on  the  specific 
control  region  mutation  motifs. 

Statistical  analysis  for  HV1  and  HVII  regions  was  made  sepa¬ 
rately,  calculating  the  nucleotide  and  sequence  diversity15  and  the 
same  parameters  were  also  calculated  for  both  regions16  using  the 
Arlequin  v.3.0  software.17  The  phylogenetic  study  was  made  by 
analysing  the  molecular  distances  between  the  studied  population 
and  chosen  populations  from  the  five  continents.18-33  Genetic 
affinities  among  populations  were  computed  using  Arlequin  v.3.0 
software  (Excoffier  et  al.,  2005).  The  phylogenetic  representation 
obtained  from  the  genetic  distances  was  made  with  the  Neighbour 
method  of  Phylip  V.3.5C  software34  and  the  tree  was  obtained  with 
Treeview  1.5.2.  software  35 


3.  Results  and  discussion 

In  this  study,  97  haplotypes  were  identified,  90  of  them  being 
unique,  expressed  by  122  polymorphic  nucleotide  positions  (Table  1 ). 
The  analysing  HV1  and  HVII  regions  separately,  67  different  haplotypes 
were  recognized  for  the  HVI  region  (55  unique)  and  71  different 
haplotypes  for  the  HVII  region  (52  unique).  The  highest  density  of 
polymorphic  sites  was  obtained  for  HVI,  which  contains  69  poly¬ 
morphic  sites  in  total  length  367  bp.  HV2  region  presents  53  poly¬ 
morphic  sites  in  268  bp.  The  most  common  polymorphism  detected  in 
the  HVI  region  was  the  single  nucleotide  substitution  T16189C,  which 
is  present  in  65.5%  of  the  sequences,  followed  by  the  T16217C,  poly¬ 
morphism  that  appeared  in  64.5%  of  the  samples.  Also  common  in  this 
Bolivian  sequences  is  the  polymorphism  in  position  16183  (62.7%)  that 
appear  as  a  transversion  A16183C  (60.9%),  or  as  a  transition  A16183G 
(0.9%),  or  deleted  (0.9%).  In  HV2  region,  beyond  the  polymorphic  site 
A263G  and  315.1C  that  appear  in  all  the  human  sequences,  the  single 
nucleotide  variant  A73G  has  a  very  high  rate  of  substitution  in  the 
sequences  (91.8%).  At  only  3  out  of  122  positions  (2.46%)  more  than 
one  mutation  events  were  noticed,  all  in  the  HVI  region.  A 16105K  (T/G) 
point  heteroplasmy  was  found  in  HVI  region  in  sequence  LPAZ038 
and  a  162Y  (T/C)  point  heteroplasmy  in  the  HVII  region  in  sequence 


Table  3 

“Estimated  haplogroup"  and  “expected  haplogroup”  designed  by  mtDNAmanager  for  sequences  LPAZ015,  LPAZ017,  LPAZ027,  LPAZ038,  LPAZ041,  LPAZ092  and  LPAZ096  based 
on  an  individual  query.  N  —  number  of  occurrence  of  the  haplotype;  Est.  Hg  —  “estimated  haplogroup”;  Exp.  Hg  —  “expected  haplogroup”. 


HVI 

HVII 

N 

Est.  Hg 

Exp.  Hg 

LPAZ015 

16168T;  16183C;  16189C;  16217C; 16270T;  16278T3 

73G;  263G;  315.1C 

1 

B4 

B4f 

LPAZ017 

16168T;  16182C;  16183C;  16189C;  16217C;  16218Ta 

73G;  263G;  315.1C 

1 

B4 

B4f 

LPAZ027 

16183C;  16189C;  16217C;  16261T;  16319Aa 

73G;  75A;  263G;  31 5.1  C;  340Ta 

1 

B4 

B4a 

LPAZ038 

16105T/G;  16168T;  16183C;  16189C;  16217C;  16301Aa 

73G;  204C;  263G;  309.1C;  315.1C3 

1 

B4 

B4f 

LPAZ041 

16066G; 16183d; 16186T;  16189C;  16217C 

73G;  143 A;  210G;  263G;  309.1C;  315.1C 

1 

H 

B4 

LPAZ092 

16181G; 16189C;  16217C 

73G;  185A;  249d;  263G;  290d;  291d; 

309.1C;  315.1C3 

1 

H 

B4 

LPAZ096 

16140C;  16183C;  16188T;  16189C;  16217C 

263G;  315. 1C 

2 

B5 

B4clb 

Means  length  heteroplasmy. 
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LPAZ061  (Table  1  ).36  Populations  with  a  high  frequency  of  unique 
haplotypes  are  consistent  with  an  early  local  settlement.5 

For  this  population  the  values  of  sequence  diversity  are  0.948 
and  0.986  for  HV1  and  HVII  regions,  respectively.  The  sequence 
diversity  when  both  hypervariable  regions  are  analysed  together  is 
0.996.  The  nucleotide  diversity  for  HVI  region  is  0.011  and  for  HVII 
region  is  0.015.  The  nucleotide  diversity  for  both  sequenced  regions 
is  0.015  ±  0.0075  (Table  2).  The  diversity  parameters  were  higher 
for  the  HVII  region  than  for  the  HVI  region  and  sequence  diversity 
increases  significantly  when  both  regions  are  analysed  together, 
which  means  that  the  analysis  of  the  HVII  region  adds  value  to  the 
population  study.  The  obtained  diversity  parameters  of  the  non¬ 
coding  hypervariable  region  I  (HVI)  are  similar  to  the  ones 
described  in  the  literature  on  Native  American  populations.37 

The  rate  of  length  heteroplasmy  is  36.36%  in  HVI  region  and 
52.73%  in  HVII  region.  It  should  be  noticed  that  the  high  rate  of 
length  heteroplasmy  that  was  found,  which  is  much  higher  than  in 
Caucasian  populations,  is  similar  to  other  populations  from 
Southwest  Asia,  namely  the  Japanese  and  Korean  populations.  The 
most  frequent  polymorphisms  found  in  this  study  are  also  frequent 
in  the  Korean  population.38,39 

It  was  possible  to  include  100%  of  the  mtDNA  sequences  into 
a  specific  mtDNA  haplogroup  according  to  nucleotide  substitutions 
of  the  HVI  and  HVII  regions.  Since  the  mtDNAmanager  is  a  dynamic 
Web-based  tool  that  follows  the  evolution  of  the  state  of  the  art,  the 
presented  results  refer  to  the  analysis  performed  on  8th  June  2009. 
The  bioinformatics  tools  of  mtDNAmanager  designate  the  “esti¬ 
mated  haplogroup”  and  the  “expected  haplogroup”.14  On  Table  1,  it 
is  presented  the  “estimated  haplogroup”  that  is  designated  when 
the  query  data  indicate  the  presence  of  accompanying  mutations 
additional  to  the  clear  diagnostic  mutations.  The  “expected  hap¬ 
logroup”  is  elected  when  the  query  data  possess  clear  diagnostic 
mutations.  The  haplogroup  distribution  (Fig.  1)  is  as  follows:  57.27% 
B4, 19.09%  Cl,  10.00%  A2,  3.64%  Dl,  2.73%  D4h3, 1.82%  H,  and  0.91% 
for  each  of  the  haplogroups  A4,  B4cla,  CZ,  D4J,  M7a  and  M8/N9b. 
The  majority  of  lineages,  96.36%,  have  polymorphisms  belonging  to 
major  haplogroups  A,  B,  C  or  D,  characteristic  of  South  American 
populations40-42  and  3.64%  of  the  lineages  belong  to  West  Eurasian 
Haplogroup  H  and  East  Asian  haplogroups  M7a  and  M8/N9b. 
Events  related  to  colonial  activities,  modern  migration  or  recent 
admixture  between  co-existing  Native  Americans  and  European 
populations  can  explain  the  appearance  of  the  haplogroup  H,  which 
is  a  minority  in  the  study.  For  sequences  LPAZ041  and  LPAZ092 
(Table  3)  the  “estimated  haplogroup"  on  mtDNAmanager  was 
haplogroup  H,  however  these  haplotypes  have  the  specific  motif  for 
haplogroup  B4  (16217C)  and  the  “expected  haplogroup”  in 
mtDNAmanager  for  these  sequences  is  haplogroup  B4.  Also,  for 
sequences  LPAZ015,  LPAZ017,  LPAZ027,  LPAZ038  and  LPAZ096 
mtDNAmanager  designated  a  different  “estimated  haplogroup”  and 
an  “expected  haplogroup”,  so  we  made  an  individual  query  for  each 
of  these  sequences  (Table  3).  The  present  haplogroup  results  are  in 
concurrence  with  previous  studies  based  on  the  mtDNA  analysis  of 
American  populations.23,37,43-46  A  study  implemented  in  native 
Bolivian  populations  showed  that  the  haplogroup  B  has  a  higher 
frequency  in  the  highlands  and  the  haplogroup  C  is  more  repre¬ 
sentative  in  the  lowlands.  In  the  Aymara  group  the  haplogroup  B 
was  the  most  frequent  and  in  Quichua  group  the  haplogroups  C  and 
A  were  more  representative.43 

By  the  analysis  of  the  phylogenetic  tree  (Fig.  2)  obtained  from 
the  molecular  distances  of  populations,  it  was  found  that  the 
Bolivian  population  reveals  lower  genetic  distance  to  the 
Korean  ($st  =  0.09976,  P  =  0.00000),  Chilean  aboriginal  Pehuenche 
($st  =  0.10402,  P  =  0.00030),  Japan  (<f>st  =  0.11130,  P  =  0.00000) 
Chilean  aboriginal  Mapuche  (<J?st  =  0.13124,  P  =  0.00000)  and 
Australian  (4>st  =  0.14335,  P  =  0.00000)  populations.  Nevertheless, 


genetic  distance  analysis  resulted  in  a  high  significant  differentia¬ 
tion  between  Bolivian  population  and  any  other  population  of  the 
comparative  study. 

Recent  studies  based  on  mtDNA  presented  evidences  for 
a  “Three-stage  model”  consisting  on  an  early  expansion  into  Beringia 
followed  by  ~  20,000  years  of  population  stability  before  the  final 
entry  into  the  Americas.47  But  in  spite  of  this  theory,  a  maritime 
entrance  in  South  America48  through  the  Pacific  Ocean  has  been 
marking  position  as  a  result  of  archaeological  findings  along  the 
different  areas  of  South  America.  On  a  recent  work,  Dillehay  et  al., 
200849  found  on  archaeological  site  Monte  Verde  in  Southern  Chile 
remains  of  nine  species  of  marine  algae  dated  between  14,220  and 
13,980  calendar  years  before  present.  Our  phylogenetic  results  — 
namely  the  lower  genetic  distance  to  Korean  —  make  us  believe  that 
the  settlement  of  South  America  was  also  along  the  Pacific  coast. 

4.  Conclusions 

Due  to  the  high  frequency  of  unique  haplotypes,  the  studied 
population  shows  a  great  genetic  variability.  The  nucleotide  and 
sequence  diversity  are  high  and  similar  to  other  populations  of 
South  America.  The  diversity  parameters  are  higher  for  HVII  region 
than  for  HVI  region  and  when  both  regions  are  analysed  together 
the  diversity  parameters  increase.  Native  Bolivians  from  La  Paz 
exhibited  haplotypes  that  belonged  to  the  haplogroups  observed  in 
Native  Americans  and  the  most  represented  is  haplogroup  B4. 
Phylogenetic  analysis  reveals  proximity  to  the  Korean,  Chilean 
aboriginal,  Japanese  and  Australian  populations,  but  the  Bolivian 
population  presents  significant  differences  when  compared  with  all 
of  the  other  chosen  populations.  Interpopulation  genetic  variability 
was  high,  suggesting  an  early  settlement  of  the  population. 
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